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Abstract 

t 2 

Some properties of the inverse of the function N(x) = -4= f_ 00 &~ 2 dt 
are studied. Its derivatives, integrals and asymptotic behavior are pre- 
sented. 

1 Introduction 

It would be difficult to overestimate the importance of the standard normal (or 
Gauss) distribution. It finds widespread application in almost every scientific 
discipline, e.g., probability theory, the theory of errors, heat conduction, biology, 
economics, physics, neural networks etc. It plays a fundamental role in 
the financial mathematics, being part of the Black-Scholes formula PJ, and its 
inverse is used in computing the implied volatility of an option Yet, little 
is known about the properties of the inverse function, e.g., series expansions, 
asymptotic behavior, integral representations. The major work done has been 
in computing fast and accurate algorithms for numerical calculations £Q. 

Over the years a few articles have appeared with analytical studies of the 
closely related error function 

X 

erf(x) = — = / e~ t2 dt 

V 71 " J 
o 

and its complement 

oo 

2 f 2 

erfc(x) = — p= / e~* dt . 
V n J 

X 

Philip ^H] introduced the notation "inverfc(x)" to denote the inverse of the 
complementary error function. He gave the first terms in the power series for 
inverfc(x), asymptotic formulas for small x in terms of continued logarithms, 
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and some expressions for the derivatives and integrals. Carlitz [3], studied the 
arithmetic properties of the coefficients in the power series of inverfc(a;). Strecok 
computed the first 200 terms in the series of inverfc(ir), and some expansions 
in series of Chebyshev polynomials. Finally, Fettis [Sj studied inverfc(a;) for 
small x, using an iterative sequence of logarithms. 

The purpose of this paper is to present some new results on the derivatives, 
integrals, and asymptotics of the inverse of the cumulative standard normal 
probability function 



which we call S(x). In section 2 we derive an ODE satisfied by S(x), and solve 
it using a power series. We introduce a family of polynomials P n related to the 
calculation of higher derivatives of S[x). In section 3 we study some properties 
of the P n , such as relations between coefficients, recurrences, and generating 
functions. We also derive a general formula for P n using the idea of "nested 
derivatives", and we compare the P n with the Hermite polynomials H n . 

In section 4 we extend the definition of the P n to n < and use them to 
calculate the integrals of S(x). We also compute the integrals of powers of 
S(x) on the interval [0,1]. Section 5 is dedicated to asymptotics of S(x) for 
x — > 0, x — > 1 using the function Lambert W. With the help of those formulas 
we derive an approximation to S(x) valid in the interval [0, 1] with error e, 
\e\ < 0.0023. Finally, appendix A contains the first 20 non-zero coefficients in 
the series of S(x), and the first 10 polynomials P n . 

2 Derivatives 

Definition 1 Let S(x) denote the inverse of 






(1) 




Proposition 2 S(x) satisfies the IVP 



S" = S(S') 2 



(2) 
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Proof. Since N(0) = 1/2, in follows that 5(1/2) = 0. From (2.1) we get 

S'[N(x)} = — = V^e^ = V^e 1 ^ 
TV (x) 

Substituting iV(x) = y we have 



Differentiating ln[6"(y)], we get (2.2). ■ 
Proposition 3 

S^ n+2 \x) =J2i2( n ) ( f )s^- i \x)S^ +1 \x)S^+ 1 \x), n > 0. 
i=o j=o ^ z ' 

Proof. Taking the n th derivative of (2.2), and using Leibnitz's Theorem, we 
have 

s( n+ v =jr( n )s( n - l Hs's')^ 

i=0 ^ l ' 

= it, (™) 5( " ~ l) XI (*) S {i ~ i+1) S U+1) 



=0 ^ ' j=0 

n i 

' n \ ( i 



■ 

Corollary 4 //£>„ = £#(±) , tfien 

£» 2n = 0, n > 0. 
Putting £>„ = (27r)^C„ , we can write 

n>0 v ; 

where 

Ci = l, C 3 = l, C 5 = 7, C 7 = 127, ... . 

Proposition 5 

5" (n) = P„-i(5)(S") n n>l 
where P n (x) is a polynomial of degree n satisfying the recurrence 

P (x) = 1, P n (x) = Pn-i{x) + nxP„_x(x), n > 1 (3) 

so i/iai 

Pi(x)=x, P 2 (x) = 1 + 2x 2 , P 3 (x) = 7x + 6x 3 , .... 
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Proof. We use induction on n. For n = 2 the result follows from (2.2). If we 
assume the result is true for n then 

S (n+1) = [Pn _^ S )(S') n ]' 

= i^_ 1 (5)5'(S') n + Pn-i(S)n(ST^S" 
= P;_ 1 (5)(5')™ +1 + P„_i(5)n(S') n_1 5(5') 2 
= [P;_ 1 (5)+n5P„_ 1 ( > S)](5'r +1 
= P„(5)(5') n+1 

Since P„_i(a;) is a polynomial of degree n — 1 by hypothesis, is clear that 

P n (x) = Pn-i(x) + nxP n --t(x) 
is a polynomial of degree n. ■ 
Corollary 6 

C„ = P„_i(0) 



3 The polynomials P n (x) 

Lemma 7 J/ we rante 

n 

Pn(x) = J2Qk xk 



we have 



Proof. 



k=0 



SO — Vl 



Q^nQri' + ^ + i)^! 1 fc=l,..-,n-2 (4) 
Qk=nQk~ 1 k = n-l,n 



Q^x fc — p„ 

fc=0 

= -j-Pn-i + nxP n -i 
ax 



x h+1 



= j2Q n k - 1 kx k - 1 + J2 n Qk- 1 

k=0 k=0 

n—1 n 

= Y,Qk+l(k + i)x k + J2 n Qk-l 



x k 



k=0 k=l 



4 



Corollary 8 In matrix form (3.1) reads A^Q"" 1 = Q n where A^ G x n 

is given by 

i, j = i + l, i = l,...,n-l 
4?i={ n > i = ^ - 1, i = 2,...,n+l 



0, 



ow 



In other words, A^ n > is a rectangular matrix with zeros in the diagonal, the 
numbers 1, 2, 3, . . . , n — 1 above if, f/ie number n below it and zeros everywhere 
else 



A (n) 



and 



10 

n 2 

n 3 















n- 1 
n 
n 



g 11 



Mtft i/ie fteip o/ these matrices, we have an expression for the coefficients 
ofP n {x) 



Q n — Y[ A ( "~ fe+1) = A (n) A (n_1 ' • • • A (1 



k=l 



Proposition 9 The polynomials P n {x) satisfy the recurrence relation 



Proof. 



— — \ 2 / V 7 

i=0 j=0 v 7 v 



/\n+2 



S (n+2) = EE ( j r.)s (n_i) s (i_j+1) s (j ' +1) 



— ' V l / V 7 



i=0 j=0 



= EE " r> I1 -,-i(S)(5T^,- J (5)(S')'- ,+1 ^(S)(S') ,+1 



i=0 j=0 



(st +2 EE ( "J (^jPn-i-mp^wPjis). 
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Proposition 10 The exponential generating function of the polynomials P n (x) 
is 

= F{x,t) = Y,Pk{x) tk 

k>Q 



}S 2 [N(x)+tN'(x)]- 



k\ 



Proof. Since 

F(M) = £ft(*£ 
fc>o K - 

k>l k>l 

J +fe+l +fe+l 

= 1 +/ 7^0^ s)ds + artF(x,i) 
Jo o s 

it follows that i) satisfies the differential-integral equation 

1 + (art - l)F(a;, t) + — F(x, s)ds = 
cte Jo 

Differentiating (3.2) with respect to t we get 

xF(x, i) + (art - 1) J^(*, t) + J^(z, t) = 



whose general solution is of the form 

F(x, t) = e'^G [te-^r + ^ 



N(x) - 



for some function G{z). 

From (3.2) we know that F(x,0) — 1, and hence 



G ( 



N{x) 



= e 2 , 



which implies that 
Therefore, 



G(z) = e^vfc + O. 
F{x,i) = e^ N ^ +tN '^-^ 



(5) 



G 



Definition 11 We define the "nested derivative" by 
3>W[/]0r) = l 



Example 12 1. 



71 > 1 



5. 



\x] = 1 



5. 



Proposition 13 



3)W [x 2 ] = (n+ 



P n (ar) = e-f^DW 



e 2 



Proof. We use induction on n. For n = the result follows from the definition 
of 2)W. Assuming the result is true for n — 1 

P n (x) = Pn-i(x) + nxP n -i[x) 

= ±\e-^* +nxe-^ llx2 & n - 1 \e^) 
dx 

= -(n - l)xe-^ x2 & n - 1 \e^)+e-^ x2 ^[V {n) {e^)] + 

dx 

+ nxe- {J! ^ llx2 & n - 1 \e^) 

c £(«-i)( e 4) + A S («-i)( e #) 
dx 



(i-l) „2 12 



d r i 



f/.T 



e 2 "* S)< n_1) (e 



e 2 



Summary 14 M^e conclude this section by comparing the properties of P n (x) 
with the well known formulas for the Hermite polynomials H n (x) ITy. Since 
the H n are deeply related with the function N(x), we would expect to see some 
similarities between the H n and the P n . 
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Pn{x) 


H n (x) 


P n (x) = e-^D( n \e^) 


H n (x) = (-l)"e x2 £;(e-* 2 ) 


J2 P k (x)£ = e hS 2 [N{x) + tN'{ x )]-*± 
k>0 


fe>0 


P n (x) = P^_i(z) +nxP n - 1 (x) 


H n (x) = -K^ix) + 2xH n _ 1 (x) 


= P n -i{S)(S') n 





4 Integrals of S(x) 



Definition 15 

S ( -~ n \x) = J j ■■■ j S(x n ) dxndxn-x . . . dxi, n > 1. 



x x± x n -i 







Lemma 16 



P n - 1 {x) = e-^ 



P„_i(0) + j e? t2 P n (t)dt (6) 
o 

Proof. It follows immediately from solving the ODE for P n -i in terms of P n . 
■ 

Proposition 17 Using (4-1 ) to define P n (x) for n < yields 



P-i(x) = x 

P-2(X) = -1 

P_ 3 (x) = -V^e x2 N (V2x) 



(7) 



and the relation 



=p n _ 1 (s)(sy 



still holds. 
Proof. For n = we have 



P_i(z) =P-±(0)+x 
S = 5 (0) = P-i(S) 
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so P-i(O) = 0. For n = -1 

P_ 2 (x) [P_ 2 (0) + 1]-1 

We can calculate S^~^ explicitly by 

x S(x) 

S(~ 1 \x) = js{t)dt = -^= J S[N{z)]e-^dz 



S{x) 



ze 2 dz 



Hence, P_ 2 (0) = -1. 
Finally, for n = — 2 



P- 3 (x) = 



P_ 3 (0) - yftN(y/2x) + ^ 



A similar calculation as the one above, making a change of variables t — N(z) 
in the integral of S < '~ 1 \x) yields 

and we conclude that P_3(0) = — m 



Corollary 18 



Proposition 19 



J S n (x)dx 
o 



n(2i + l), n = 2fc, k>l 

i=l 

0, n = 2fc + l, fc>0 



Proof. 



1 oo 

y s*™(x)dx = y z"iv'(z)d2 



z"e 2 Z dz 



Y[(2i + 1), n = 2fc, fc>l 



i=l 
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5 Asymptotics 

Definition 20 We'll denote by LW(x) the function Lambert W 

LW{x)e LW{x) = x 
This function has the series representation 

LW(x) = £ M^V, 

n>l 

the derivative 

and it has the asymptotic behavior 

LW(x) ~ ln(x) — ln[ln(x)] x — > oo. 

Proposition 21 



5(1) ~ g (x) = -WW ) ' •'• - !! 



S(x) ~ gi{x) = x j LW ( ), x^l 



1 

2n(x - l) 1 



Both functions go(x) and <?i(x) satisfy the ODE 

2 



t" = s(sT 



1 



3 2 (l + .9 2 ) 



fffoT, lor \g\ 



Proof. 



1 1 

N(x) ~ e 2 — , x — * — oo 

V27T X 

1 S(t) 2 1 

V21F 5(*)' 

, , s<t) 2 1 

5te— =-, i^O 



Using the definition of LW(x) we have 
or 



The case x — > 1 is completely analogous. 
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Corollary 22 Combining the above expressions, we can get the approximation 

S(x) * g 2 (x) = (2x - 1)^( 2 ^_ 1)2 ) W 

good through the interval (0, 1). 

VFe can re/me it even more by putting 

S(x) * 33 (x) = Q {X ^L W (^1-^) (10) 
Q(x) = -1 + (6 - 2\/2~T)a; + (-12 + 6V2^)x 2 + (8 - 4V2~r)a; 3 
w/iere Q(x) /ias been chosen such that 

Q(0) = -1, Q(l) = 0, Q(l/2) = 0, Q'(l/2) = v^, Q"(l/2) = 

6 Appendix 

The first 10 P„(x) are 

P (x) = 1 
Pi(x) = x 
P 2 ix) = l + 2x 2 
P 3 (x) = 7a; + 6x 3 
P 4 {x) = 7 + 46a; 2 + 24a; 4 
P 5 (x) = 127a; + 326a; 3 + 120a; 5 
P 6 (x) = 127 + 1740a; 2 + 2556a; 4 + 720a; 6 
P 7 (x) = 4369a; + 22404a; 3 + 22212a; 5 + 5040a; 7 
P s (x) = 4369 + 102164a; 2 + 290292a; 4 + 212976a; 6 + 40320a; 8 
P 9 (x) = 243649a; + 2080644a; 3 + 3890484a; 5 + 2239344a; 7 + 362880a; 9 
P w (x) = 243649 + 8678422a; 2 + 40258860a; 4 + 54580248x 6 + 
+ 25659360a; 8 + 3628800a; 10 

The first few odd C n are 
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n 



1 

3 

5 

7 

9 

11 

13 

15 

17 

19 

21 

23 

25 

27 

29 

31 

33 

35 



108026349476762041127839800617281 
76683701969726780307420968904733441 



15773461423793767 
4591227123230945407 
1598351733247609852849 
655782249799531714375489 
313160404864973852338669783 
172201668512657346455126457343 



243649 
20036983 
2280356863 
343141433761 
65967241200001 



1 
1 
7 

127 

4369 



37 61154674195324330125295778531172438727 
39 54441029530574028687402753586278549396607 
41 53789884101606550209324949796685518122943569 
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